Superconductivity in the Orbital Degenerate Model for Heavy Fermion Systems 
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Magnetism and superconductivity of new heavy fermion compounds CeTIns (T=Co, Rh and Ir) 
are investigated by applying fluctuation exchange approximation to an orbital degenerate Hubbard 
model. The superconducting phase with d^2_j^2 -symmetry is found to appear next to the antiferro- 
magnetic phase with increasing the orbital splitting energy. The present theory suggests that the 
orbital splitting energy plays a key role of controlling parameter for the quantum phase transitions 
in the heavy fermion system. 
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I. INTRODUCTION 

Superconductivity in strongly correlated electron sys- 
tems has been one of central issues in the research field 
of condensed matter physics, since the pioneering dis- 
covery of superconductivity in CeCu2Si2 0]. The sub- 
sequent discovery of high temperature superconductivity 
in cuprates has accelerated further investigations of this 
subject, leading to unambiguous identification of the un- 
conventional nature of dx2_y2--wave pairing in cuprates. 
Recently, new heavy fermion compounds CeTIns (T=Rh, 
Ir, and Co) have been discovered ||]. Among them, 
CeRhlus exhibits an antiferromagnetic (AFM) transition 
at a Neel temperature Tn^S.SK and becomes supercon- 
ducting only under hydrostatic pressure larger than 15 
kbar. On the other hand, both Celrlns and CeCoIns 
are superconducting at ambient pressure with transition 
temperatures rc=0.4K and 2.3K, respectively. Espe- 
cially, the Tc of CeCoIns is the highest among those yet 
observed for heavy fermion superconductors H . 

Reflecting the fact that CeTIns has the HoCoGas-type 
tetragonal crystal structure, quasi two-dimensional Fermi 
surfaces have been observed in de Haas-van Alphen ex- 
periments for the compounds, consistent with the band- 
structure calculation results [Q. Concerning the super- 
conducting state, nuclear relaxation rate in CeTIns ex- 
hibits behavior below Tc ^ and thermal conductiv- 
ity in CeCoIns is found to include a component with 
four-fold symmetry strongly suggesting d^2_j^2-wave 
pairing in the superconducting phase of CeTIns. Fur- 
thermore, it has been shown that in the alloy system 
CeRhi_2:Ira;Ins, the superconducting phase appears in 
the neighborhood of the AFM phase . A natural con- 
sequence of these experimental results is that supercon- 
ductivity in CeTIns compounds is induced by AFM spin 
fluctuations, similarly to high-Tc superconductivity. 

In spite of the similarities mentioned above, here 
we should emphasize several differences between heavy 
fermion superconductors such as CeTIns and high-Tc 
cuprates. First of all, relevant electrons in high-Tc su- 
perconductors are almost itinerant 3d-electrons, while 
for heavy fermion superconductors, relevant ones are /- 
electrons, which almost localized and their dispersion 



is mainly determined by hybridization with conduction 
electrons. Another important difference is concerned 
with the electronic states relevant to low-energy physics. 
For cuprates, it is widely recognized that a single-band 
model for holes in d^^-y-^ orbitals on the square lattice 
is a good starting point. On the other hand, for heavy 
fermion superconductors, several Fermi surfaces are ob- 
served in general. Such complex electronic states may 
be eventually traced back to the orbital degeneracy and 
relatively weak crystalline electric fleld (CEF) effect of 
/-electrons compared with the 3d-electrons in the Cu02 
plane. It implies that construction of a realistic micro- 
scopic model is not an easy work for heavy fermion sys- 
tems. Therefore, up to now, most of theoretical investi- 
gations for superconductivity in heavy fermion systems 
have been restricted in a phenomenological level. 

In this paper, we discuss the effect of orbital fluctua- 
tions on superconductivity based on a microscopic theory 
applied to a microscopic model. In the next section, we 
introduce an orbital degenerate model obtained including 
important characters of CeTIns . Then, in order to study 
the superconducting transition induced by spin and/or 
orbital fluctuations, we develop a strong-coupling the- 
ory using fluctuation exchange (FLEX) approximation 
in which spin and orbital fluctuations as well as the 
single-particle spectrum are determined self-consistently. 
Finally, we discuss experimental results for CeTIns in the 
light of the present theory. 



II. MODEL HAMILTONIAN 

In order to introduce a minimal model for /-electron 
systems, we start from the local basis for a Ce'^^ ion. 
Among 14-fold degenerate 4/-electronic states, due to 
the effect of strong spin-orbit coupling, only the j=5/2 
sextet effectively contributes to the low-energy excita- 
tions (j is total angular momentum). Under the cu- 
bic CEF, this sextet is further split into Fy doublet 
and Fg quartet, where the corresponding eigen-states are 
given by |r7±)= ^\ ± 5/2)-^\ T 3/2), \Ti2) = 

y/5/6\ ± 5/2) + ^lyel T 3/2), and |r^^2)=l ± 1/2)- Here 
+ and — in the subscripts denote "pseudo-spin" up and 
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down in each Kramers doublet, respectively. For sim- 
plicity, we consider only \Tg,) states, further split into 
two Kramers doublets as IFg^) and \T^2) under the addi- 
tional tetragonal CEP. Note that |r^2) ^-nd |rf2) belong 
to Tj and Fg irreducible representations, respectively, in 
the tetragonal system. Validity of this assumption for 
CeTIn5 will be discussed later in this paper, but here we 
stress that the Hamiltonian constructed from Fg quar- 
tet is the simplest model including essential physics of 
interplay between spin and orbital fluctuations. 

In order to consider itinerant features of 4/-electrons, 
we take into account nearest-neighbor hopping of /- 
electrons by the tight-binding method It should be 
noted that the matrix elements of the nearest-neighbor 
hopping depend on not only the orbital but also the hop- 
ping direction, since the forms of wavefunctions of |Fg2) 

(2) 

and iFgji) states are different from each other. We can 
determine the hopping matrix elements by assuming that 
the overlap integral through the cr-bond (ffcr) is dom- 
inant. Since CeTIns has a tetragonal crystal structure 
and quasi two-dimensional Fermi surfaces have been ex- 
perimentally observed |Q,|lO| , it is natural to consider the 
two-dimensional square lattice composed of Ce'^"'" ions. 
Although the hybridization of / electrons with In 5p elec- 
tronic states may be the main source of kinetic energy for 
/ electrons, in the present scheme such an effect is con- 
sidered as the effective hopping amplitude for / quasi- 
particles, after the p-electron degrees of freedom are in- 
tegrated out. 

By further adding the on-site Coulomb interaction 
terms among /-electrons, the effective Hamiltonian with 
orbital degeneracy is given by 

H= 0/L/i+arV-eE("il--"i2.)/2 
iarr'cj i 

+ U'Ynirinirl+U'Yniicrni2a', (1) 
ir iaa' 

where /iro- is the annihilation operator for an /-electron 
with pseudo-spin a in the r-orbital state Fg at site i, 
a is the vector connecting nearest-neighbor sites, and 
"'ir(T=//T-o-/iTCT- The first term represents the nearest- 
neighbor hopping of /-electrons with the amplitude t^^, 
between r and r' orbitals along the a-direction, given by 
q^=-^/3q^= -V3tfi=3tj2=l for a=x and tli=y^tl2 = 
\/3i2i— 3^22=1 for ^—Y: respectively, in energy units 
where tii—1. The second term denotes the tetragonal 
CEF, leading to an energy splitting e between the two 
orbitals. In the third and fourth terms, U and U' are 
the intra- and inter-orbital Coulomb interactions, respec- 
tively. Due to the rotational invariance in the orbital 
space for the interaction part, U' should be equal to U, 
when we ignore the Hund's rule coupling, since it is ir- 
relevant in the large-?/ limit for the quarter-filling case 
with one /-electron per site. Thus, in this paper, we re- 
strict ourselves to the case of U=U'. Note also that in 



the quarter-filling case, the present model is virtually re- 
duced to the half-filled single-orbital Hubbard model in 
the limit of e=oo. 



III. FLEX APPROXIMATION 

In our previous work, we have developed a weak- 
coupling theory for superconductivity based on the same 
orbital degenerate model described above, using the 
static spin and orbital fluctuations obtained within the 
random phase approximation (RPA) Various su- 

perconducting phase have been found around varieties of 
ordered phases, whose boundaries are determined by the 
RPA instability. In order to develop more sophisticated 
theory, we should include effects of (1) mode- mode cou- 
pling and (2) quasi-particle damping, neglected in the 
previous work. Regarding (1), the RPA does not incor- 
porate effects of mode-mode coupling between fluctua- 
tions, but the mode- mode coupling modifies significantly 
the temperature and frequency dependences of spin and 
orbital fluctuations. Concerning (2), within the weak- 
coupling theory for superconductivity, damping of quasi- 
particles by the scattering due to spin and orbital fluc- 
tuations is not taken into account, which leads to the 
suppression of superconductivity. Therefore, the weak- 
coupling theory generally overestimates the region of su- 
perconducting phase. 

In the present paper, we apply the fluctuation ex- 
change (FLEX) approximation ||] to the orbital degener- 
ate model discussed in the preceding section. The FLEX 
approximation has the following two features: (1) This 
is a kind of the mode-mode coupling theory, where spin 
and orbital fluctuations and the spectra of /-electrons 
are self-consistently determined through the fluctuation 
exchange self-energy. (2) It provides the Dyson-Gorkov 
equation where the normal and anomalous self-energies 
are obtained on an equal footing. Here we emphasize that 
the FLEX approximation has been successful to under- 
stand consistently the normal and superconducting states 
of high- Tc cuprates, in particular with reasonable estima- 
tion of Tc |12|. Therefore, it is interesting to apply the 
FLEX approximation to the orbital degenerate model to 
understand its properties concerning superconductivity. 

In the doubly degenerate case, the Green's functions 
for /-electrons form a 2x2 matrix and they follow the 
Dyson-Gorkov equations 



(2) 



where G*^"-' (k) is a matrix for the non-interacting Green's 
function. Within the FLEX approximation, components 
of the self-energy matrix (k) are given by 



^^lik) = ^ E E Kl^ii<l)G,Ak - q), 



(3) 



with 
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+ U°)m{W + U°) + - lu\,nM, (4) 

where the first and second terms of (q) give contri- 

butions of the spin and orbital fluctuations, respectively, 
to the self-energy. 

In Eq. (4), )<f{q) and x°{q) a^re the 4x4 matrices of the 
spin and orbital fluctuations, given by 

r{q) = [i - U'f{q)]-'il{q), (5) 

r{q) = [i+u°mr'f{q), (6) 

where the matrix element of the irreducible susceptibility 
Xij,st{q) is given by 

Xij,sfiq) = -j^'^Gstik + q)Gjtik), 



(7) 



For the susceptibility matrices, the labels of row and col- 
umn appear in the order 11, 22, 12, and 21, pairs of 
orbital indices 1 and 2. In these expressions, T is a 
temperature, Nq is number of unit cells, abbreviations 
k—{k,iu:n) and (7=(q, if2„) {cjn={2n+l)nT, rj„=2n7rT) 
are used. The interaction matrices and U° are given 
by 
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By solving the self-consistent equations, the spin and 
orbital fluctuations and the Green's functions for /- 
electrons are determined simultaneously. 

In order to discuss superconductivity, it is necessary 
to calculate anomalous self-energy. The matrix ele- 
ment of the anomalous self-energy (fc) is obtained 
by the functional derivative of the thermodynamic po- 
tential with respect to the anomalous Green's function 
as follows 



T 

Wo 



q), 



with 



F{k) = G(fc)E(2)(fc)G'(-fc), 



(8) 



(9) 



where the matrix elements of the effective pairing interac- 
tions for spin-singlet and spin-triplet channels are given, 
respectively, by 



V^^,iAq) = [-\u^r{q)U'' + \u°r{q)U° 



Vjrn.lAq) 



[lu'x\q)U' + lu°r{q)U° 



fi7n,liy ■ 



(10) 



(11) 



The Tc is obtained by the temperature at which the max- 
imum eigenvalue of Eq. (8) becomes unity. As already 
pointed out in the previous works based on the RPA 
[|l3| , p| , one can see from these effective pairing inter- 
actions that developments of both spin and orbital fluc- 
tuations have a destructive interference for the singlet 
channel, while they are constructive for the triplet one. 



IV. CALCULATED RESULTS 

The FLEX calculation is numerically carried out for 
each value of e at fixed parameter values of t/=C/'=4 and 
n—1 corresponding to one /-electron density per site. All 
summations involved in the above self-consistent equa- 
tions are performed using the fast Fourier transformation 
algorithm for the k-space with 32 x 32 meshes in the first 
Brillouin zone and for Matsubara frequency sum with 
energy cut-off five times larger than the relevant band 
width. In Fig. 1, q dependences of the principal com- 
ponents of x^(q, 0) and x°(q, 0) are shown for a fixed 
temperature T—Q.02 for different level splitting, e: The 
upper panel for e=0 and the lower panel for s—2. For 
e=0, corresponding to the orbitally degenerate case, the 
AFM spin fluctuation in the r=l orbital is enhanced, but 
not sufficiently developed to induce d^2_y2-waMe super- 
conductivity. With increasing the orbital splitting energy 
to e=2, the AFM spin fluctuation for the t=1 orbital 
further develops, and orbital fluctuations are completely 
suppressed, in comparison with the developed AFM spin 
fluctuation. 
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FIG. 1. (a) Spin and (b) orbital susceptibilities in q space 
for e=0. (c) and (d) are for e=2. 

In Fig. 2, the phase diagram obtained within the 
FLEX approximation for the orbital degenerate system 
are shown, where the solid and open circles describe 
the superconducting and AFM transitions, respectively. 
The dotted curve is a schematic phase boundary ex- 
pected between the two phases. From Fig. 2, we can 
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see that (1) the spin-singlet superconducting phase with 
Big-symmetry appears next to the AFM phase and (2) Tc 
is enhanced with increasing the orbital splitting energy 
e. From these observations, we can conclude that the 
superconducting phase is induced by the development of 
the AFM spin fluctuations for the r=l orbital with in- 
creasing the orbital splitting energy e. 
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FIG. 2. Phase diagram in the T-e plane for (7=I7'=4.0 ob- 
tained by the FLEX approximation. 



will be an interesting future problem to elucidate the role 
of orbital splitting for the more realistic two r7 model, 
which includes not only the nearest neighbor hopping 
but also the next nearest neighbor one and thus, may 
reproduce the realistic electronic states better than the 
Fg model 

In order to construct a whole story for CeTIns com- 
pounds, it may be required to use the even more realistic 
f-p model including the f-p hybridization explicitly. In 
the f-p model, we can have two different energy scales: 
The large energy scale is the band width for the conduc- 
tion electron and/or the Coulomb interaction, while the 
small energy scale corresponds to the energy splitting of 
the CEF levels discussed here. The present study indi- 
cates the possibility that the small energy scale of the 
CEF level splitting plays a key role as the controlling 
parameter of quantum phase transitions. 

In summary, based on the effective microscopic model 
with orbital degeneracy for /-electron systems, we have 
proposed that the orbital splitting energy is the parame- 
ter controlling the change from the paramagnetic to the 
AFM phase with the dx2_y2-wave superconducting phase 
in between. Actually, the self-consistent FLEX approx- 
imation applied to the orbital degenerate model shows 
that the dx2_y2-wa.ve superconducting phase is induced 
by increasing the orbital splitting energy. 



V. SUMMARY AND DISCUSSION 

Let us now discuss the experimental results for CeTIns 
in comparison with the present theoretical results. First 
we should note that in the actual tetragonal crystal, the 
four-fold degenerate Fg states in the cubic notation split 
into two Kramers doublets, F7 and Fg, as we mentioned 
in the Sec. 2. Therefore, in the tetragonal system, the 
j=5/2 states split into two F7 and one Fg. Analyses of 
experimental data of magnetic susceptibility of CeTIns 
by using the CEF theory seem to be consistent more 
with the level scheme where two F7 are lower than the 
Fg [0. The energy splitting between the two F7 is esti- 
mated as 68K for CeRhlns, 61K for Celrlns, and 151K 
for CeCoIus. According to the present analysis, higher 
Tc is obtained for larger e, consistent with the tendency 
in Tc of the two superconducting materials: rc=0.4K for 
Celrlns and Tc=2.3K for CeCoIng. 

However, the antiferromagnetically ordered CeRhIn5 
with Tn=3.8K has an intermediate value for e. One pos- 
sible scenario to understand the discrepancy is to con- 
sider difference of quasi-two dimensionality, as we have 
pointed out in To make more direct and quantita- 
tive comparison with experimental results, however, one 
should be aware of an assumption for the model used 
in the present study, which may be called as Fg model. 
Namely, this Fg model assumes that F7 and Fg are the 
lower two Kramers doublets, which may be different from 
the level scheme obtained from the experiments. Thus, it 



ACKNOWLEDGEMENT 

The authors would like to thank T. Maehira for many 
helpful discussions. T.H. and K.U. are supported by the 
Grant-in- Aid for Scientific Research from Japan Society 
for the Promotion of Science. 



REFERENCES 



[1] Steglich F et al. 1979 Phys. Rev. Lett. 43 1892 
[2] Hegger H et al. 2000 Phys. Rev. Lett. 84 4986 
Petrovic C et al. 2001 Europhys. Lett. 53 354 
Petrovic C et al. 2001 J. Phys.: Condens. Matter 13 L337 
[3] Recently it has been discovered that PuCoGas exhibits 
superconductivity with T'c = 18.5K. See Sarrao J h et al. 
2002 Nature(London) 420 297 
[4] Haga Y et al. 2001 Phys. Rev. B 63 R060503 
[5] Kohori Y et al. 2000 Eur. Phys. J. B 18 601 
Zheng G-q et al. 2001 Phys. Rev. Lett. 86 4664 
[6] Izawa K et al. 2001 Phys. Rev. Lett. 87 057002 
[7] Pagliuso P G e£ al. 2001 Phys. Rev. B 64 R100503 
[8] Bickers N E et al. 1989 Phys. Rev. Lett. 62 961 



4 



[9] Maehira T et al. 2002 J. Phy s Soc. Jpn. Suppl A 71 285 
and see also Maehira T et al. cond-mat/021136C 
Settai R et al. 2001 J. Phys.: Condens. Matter. 13 L627 
Takimoto T et al. 2002 J. Phys.: Condens. Matter 14 L369 



[10] 

[11] 
[12] 



Scalapino D J 1995 Phys. Rep. 250 329 
Moriya T and Ueda K 2000 Adv. Phys. 49 555 
[13] Takimoto T 2000 Phys. Rev. B 62 R14641 
[14] Takeuchi T et al. 2001 J. Phys. Soc. Jpn. 70 877 

Shishido H et al. 2002 J. Phy s. Soc. Jpn. 71 162 
[15] Hotta T and Ueda K 2002 |cond-mat/0211467 



5 



